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Abstract. The results of a comprehensive investigation into the effects of toroidicity- 
induced Alfven eigenmodes (TAE) and energetic particle modes on the NBI-generated 
fast-ion population in MAST plasmas are reported. Fast-ion redistribution due 
to frequency-chirping TAE in the range 50 kHz to 100 kHz, and frequency-chirping 
energetic particle modes known as fishbones in the range 20 kHz to 50 kHz, is observed. 
TAE and fishbones are also observed to cause losses of fast ions from the plasma. 
The spatial and temporal evolution of the fast-ion distribution is determined using 
a fission chamber, a radially-scanning collimated neutron flux monitor, a fast-ion 
deuterium alpha spectrometer and a charged fusion product detector. Modelling 
using the global transport analysis code Transp, with ad hoc anomalous diffusion 
and fishbone loss models introduced, reproduces the coarsest features of the affected 
fast-ion distribution in the presence of energetic-particle-driven modes. The spectrally 
and spatially resolved measurements show however that these models do not fully 
capture the effects of chirping modes on the fast-ion distribution. 
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The Mega-Amp Spherical Tokamak (MAST) was reliant on neutral beam heating 
to access high-performance operating regimes m ei- One of the world’s largest 
spherical tokamaks to date, typical parameters of MAST were: major and minor radius 
R = 0.95 m, a = 0.60 m; plasma current / p = 400 — 900 kA; toroidal field on axis 
.Bt — 0.40 — 0.58 T; core electron density and temperature n e o = 3 x 10 19 m -3 and 
T e o = 1 keV. The neutral beam injection (NBI) system installed on MAST consisted of 
two positive ion sources capable of accelerating deuterons to energies of 75keV. Each 
injector could deliver up to 2.5 MW of power when tuned to optimum perveance [3]. 

Fast ions (FI) in MAST were produced by ionisation of beam neutrals as they 
propagated through the plasma. The relatively weak magnetic field resulted in large 
Larmor radii of up to 10 cm for the highest energy beam ions. Once they are deposited, 
FI interact with the background plasma, slowing down by electron drag and scattering 
off thermal ions. Ultimately, the ions may either: thermalise, reaching equilibrium with 
the background plasma; undergo fusion reactions, producing a triton and proton, or 3 He 
nucleus and neutron; become re-neutralised by charge exchange with beam or thermal 
neutrals; collisionally scatter onto unconfined orbits; or be lost from the plasma due 
to instabilities or static field perturbations. In order to maximise the electron and ion 
temperature and fraction of NBI-driven current in the plasma, and hence measures of 
performance such as plasma beta, it is desirable to confine the FI for as long as possible. 
Maximising the fraction of NBI-driven current requires these ions to be confined at 
low plasma density, posing challenges for avoiding energetic particle-driven instabilities. 
Diagnosing the behaviour of FI in MAST and establishing the mechanisms for their 
redistribution within and loss from the plasma are the objectives of the present study. 

This paper focuses specifically on transport of FI in MAST caused by MHD 
instabilities which resonate with them [3]. This resonant interaction allows significant 
energy transfer from the FI population to the MHD mode in the linear growth phase. 
Subsequent nonlinear evolution of the mode is typically associated with large-scale 
coherent transport of the resonant FI [5]J. The mode acts to flatten the FI distribution 
close to the resonance in real space and in velocity space. In the worst-case scenario, 
where several resonances overlap or where sufficient free energy is present in the 
FI distribution to drive the mode over a spatially-extended region, FI may even be 
transported to the plasma boundary and lost. Such losses, as well as degrading plasma 
performance, may constitute unacceptable localised heat loads on the plasma facing 
components in a fusion reactor [6j. 

MAST provided a useful facility for the study of energetic-particle physics for a 
number of reasons. Foremost among these were the super-Alfvenic beam ions, which 
allowed Alfven eigenmodes to be driven via the fundamental resonance v\\ = va, and the 
large orbit widths of the FI which pose a challenge to the verification of numerical models 
[7]. MAST was therefore equipped with a number of advanced diagnostics designed to 
probe the dynamics of the beam-ion population. Recent work by many of the present 
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authors |Sj has highlighted the complementary nature of the different FI diagnostics 
installed on MAST, presenting an overview of the observed FI transport in the presence 
of chirping toroidicity-induced Alfven eigenmodes (TAE) and fishbones, sawteeth and 
the long-lived internal kink mode. This paper forms a natural continuation of that study, 
and aims to provide a quantitative verification of the key results presented in the prior 
work. The main findings are the confirmation that anomalous transport and losses 
of FI accompany both chirping TAE and fishbones; the demonstration that multiple 
diagnostic signals may in some cases be modelled consistently using anomalous FI 
diffusion in a global transport model; and the identification of the passing FI population 
in the core of the plasma as being susceptible to redistribution in the presence of 
fishbones. This work expands upon the previous study by providing: 

• a rigorous quantitative assessment of the degree to which the reduction in FI 
confinement is correlated with TAE and fishbones; 

• direct experimental confirmation of the fact that these MHD modes cause losses of 
FI from the plasma; 

• an improved implementation of the fishbone loss model in Transp, including fast- 
ion deuterium alpha (FID A) data as a constraint for the first time. 

Section [2] of this paper describes the FI diagnostics used in this study. The nature of 
the energetic-particle-driven MHD activity commonly observed in MAST is illustrated 
in section [3} Section [4] outlines the plasma scenarios investigated in the course of this 
study, and sections [5] and [6] present the experimental results and associated FI transport 
modelling. This modelling was performed with the 1.5-dimensional (ID profiles mapped 
onto a 2D equilibrium grid), time-dependent transport analysis code Transp [5j. The 
Monte Carlo particle-tracking module Nubeam m was used to model the deposition 
and orbits of FI and to generate local FI distributions as a function of energy and pitch. 
The main conclusions of the work are presented in section [7| 

2. Fast-ion diagnostics on MAST 

2.1. Description of the diagnostics 

During the 2013 MAST experimental campaign, four FI diagnostics were available on the 
device; one of these, the charged fusion product detector (CFPD), is a novel prototype 
diagnostic which was only installed for a limited period m- The other three diagnostics 
are: a uranium-235 fission chamber for volume-integrated neutron flux measurements 
U23; a collimated, radially-scanning neutron detector (known as the neutron camera, or 
NC) [13]; and a dual-view Fast-Ion Deuterium Alpha (FIDA) diagnostic [H]. Figure 
[T] shows the detector positions and lines of sight of each of these diagnostics projected 
onto the poloidal and equatorial planes of the MAST vessel. 

It should be noted that although measurements of the neutron flux are 
not fundamentally restricted to FI diagnostic applications, the relatively low ion 
temperature in MAST plasmas (typically 7] < 2 keV in the core of the plasma) resulted 
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Figure 1: Poloidal (a) and equatorial (b) projections of the viewing geometry of each 
of the MAST FI diagnostics. Flux surfaces are superimposed as dashed lines; solid 
grey and black components respectively are cross-sections of poloidal field coils and 
graphite armour on the divertor and beam dumps. Green and grey shaded areas 
show the regions spanned by toroidally- and vertically-oriented SS beam-viewing FIDA 
lines of sight (LOS) respectively. Red and blue lines show examples of NC LOS, 
both midplane and diagonally-oriented chords. Coloured, curved lines show fusion 
proton/triton trajectories from the plasma which reach the CFPD. NBI footprints are 
shown as red shaded regions. The cyan shaded area marks the region spanned by 
toroidal FIDA reference views, used for background subtraction. Also shown in (b) are 
the positions of the fission chamber (FC) and an outboard midplane Mirnov pick-up 
coil, used to detect MHD activity in the plasma. Midplane intersection radii of the 
active (SS beam-viewing) and reference (background subtraction) 

FIDA chords are shown as dark brown lines in (b). Adapted from 
3335/57/1/014006. ©IOP Publishing. Reproduced by permission of IOP Publishing. 
All rights reserved. 


vertically-oriented 
doi:10.1088/0741- 


in a very low contribution of thermonuclear fusion reactions to the total neutron yield. 
Approximately 98% of D-D fusion reactions in MAST occurred either between a beam 
ion and a thermal ion ©80%), or between two beam ions ©20%). As a result, the flux 
of 2.45 MeV neutrons leaving the plasma provided information about the non-thermal, 
high-energy part of the deuterium ion velocity distribution. The 235 U fission chamber, or 
FC H2], measured the volume-integrated neutron flux from the plasma with a temporal 
resolution of 10 ps, while the NC had four collimated lines of sight (LOS) allowing 
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the detectors to view the plasma through a thin flange in the MAST vacuum vessel 
m- Two of the LOS lay in the machine midplane, and the other two were oriented 
diagonally downward such that at the point of their tangency to the flux surfaces they lay 
approximately 20 cm below the midplane. The NC assembly was mounted on a curved 
rail, allowing the LOS to be scanned across the plasma between discharges; repeated 
discharges allowed a neutron emission profile, integrated within the NC fields of view, 
to be obtained as a function of viewing tangency radius. Data are typically integrated 
over 1 ms, setting the temporal resolution of the measurements. 

The CFPD was installed on MAST to provide a complementary measurement to 
that provided by the NC nu. Four detectors measured the flux of 3.02 MeV protons 
and 1.01 MeV tritons (the gyro-radii of these two fusion products being very nearly 
equal) produced by D-D fusion reactions. Collimators in the protective casing allowed 
fusion products to reach the detectors, and an aluminium foil prevented detection of soft 
X-rays. The diagnostic was mounted on the end of a linear reciprocating probe inserted 
into the MAST vessel at the midplane. This diagnostic technique is practical in small, 
low-field devices such as MAST because the large gyro-radii of the fusion products render 
them essentially unconfined by the magnetic field; they leave the plasma within a single 
gyro-orbit. With a suitable equilibrium reconstruction, orbits of these particles may be 
tracked ‘backward in time’ from the detector to determine the possible location of their 
emission. Retraction or insertion of the probe between repeated discharges allowed a 
spatial scan of the proton/triton emission profile. Protons and tritons are discriminated 
in post-processing. As in the case of the NC, the integration time of the data is typically 
set to 1 ms. 

The FIDA diagnostic provided a measurement of the FI density by observing the 
Balmer-alpha (D a ) emission of reneutralised fast deuterons which had undergone charge 
exchange (CX) with beam or halo neutrals [Hj. Halo neutrals are thermal neutrals 
created by CX between beam neutrals and thermal ions. The light was coupled via 
optical fibres to a spectrometer, and spectral resolution was provided by dispersion 
of a bandpass-filtered portion of this light. Spatial resolution was obtained by using 
individual fibres for each LOS through the plasma. The spectrum of FIDA light provides 
information about the distribution of the LOS velocity of the FI at the point at which 
CX occurred; light from reneutralised ions moving away from the lens is redshifted, 
while that from ions moving toward the lens is blueshifted. Eleven channels could be 
connected to the spectrometer at any given time, allowing a radial profile of the FIDA 
emission to be acquired within a single discharge. Data were acquired with 0.3 ms 
temporal resolution. Two sets of views were available: a near-toroidal set from a lens 
mounted just above the vessel midplane, and a near-vertical set from a lens mounted 
inside the vessel, looking vertically downward at the neutral beam. The toroidal views 
were predominantly sensitive to passing FI, with a large pitch p = v\\/v, while the 
vertical views were more sensitive to trapped FI with smaller values of pitch. 

An example of the data from each of the diagnostics described in this section is 
shown in figure|2| 2.0 MW of neutral beam heating power (at E \>i = 70 keV) was applied 
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from 0.051 s, with an additional 0.7 MW (at E h2 = 44keV) applied from 0.101 s. A series 
of large fishbones occurred between 0.20 s — 0.25 s; the time of the first of these bursts 
of MHD activity is indicated by the dashed vertical line. The effect of the fishbones is 
apparent as a significant drop in the magnitude of the signal in each time trace; this 
effect is particularly marked in the NC, FID A and CFPD data, where drops in excess 
of 25% of the pre-fishbone signal are observed. 

2.2. Velocity-space sensitivity 

Each of the diagnostics has a sensitivity to FI which varies throughout velocity space as 
well as real space. This is most apparent in the case of the fusion product diagnostics, 
since the D-D fusion cross-section is a strong function of reactant energy. Figure [3] 
shows the cross-section of the two branches of the D-D fusion reaction as a function 
of centre-of-mass energy from 5keV to 100 keV m- It is clear that the sensitivity of 
fusion product measurements is strongly weighted to the high-energy FI, especially given 
that the contribution of reactions at each energy to the total reactivity is proportional 
to av rather than simply to a. The weighting of the measurements at each point in 
velocity space is determined by averaging the reactivity over the velocity distribution of 
the target ions. This is a thermal Maxwellian distribution in the case of beam-thermal 
reactions, and a numerically-calculated velocity distribution of beam ions in the case of 
beam-beam reactions. 

An illustrative example of the velocity-space sensitivity of the FIDA diagnostic is 
shown in figure [4| FIDA measurements are sensitive only to FI in a bounded region of 
velocity space; this region is shown in panels (b) and (d) of figure [4] for two toroidally- 
viewing chords with different beam intersection radii. Variation of the sensitivity within 
this bounded region arises due to geometric factors, which determine the proportion of 
the gyro-orbit in which the LOS velocity of each FI renders it detectable in a particular 
viewing chord at a particular wavelength. Atomic physics factors also play a part, by 
determining the probability that interaction with beam or halo neutrals will give rise to 
the emission of a D a photon. This geometric part of this sensitivity function, known as 
a weight function, is determined by the geometry of the LOS, the magnetic field and the 
neutral beam. The CX probability depends on the electron temperature and density via 
their effect on the beam and halo neutral density and distribution of excited states. NBI 
energy and full, half and third-energy current fractions also affect the CX probability 
via their effect on the relative velocity between FI and beam neutrals. Observing FIDA 
light at a given wavelength corresponds to observing emission which is a product of 
the wavelength-dependent weight function and the FI distribution at the location of 
beam intersection. Stark splitting of the Balmer-alpha radiation also plays a part in 
determining the precise details of the weight function. Given the complexity of their 
construction, these weight functions must be calculated numerically in a similar manner 
to the fusion neutron or proton emissivity profiles. The general method for quantitative 
calculation of FIDA weight functions is detailed by Salewski et al. in ref. m, while the 
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Figure 2: Time traces from single channels of the four MAST FI diagnostics, from 
shot #29980. The first four panels show, from top to bottom: fission chamber count 
rate; neutron camera count rate from a midplane chord with tangency radius p\ FIDA 
spectral radiance at wavelength A from a chord intersecting the axis of the SS neutral 
beam at major radius R\ CFPD proton count rate from a channel with fusion proton 
trajectories intersecting the midplane at major radius R (where the values of p, R and 
A are given by the labels in each panel). The thickness of the lines shows the statistical 
uncertainty in the data. Shown in the bottom panel is the total NBI power and the 
RMS voltage in vertically-oriented Mirnov coil on the outboard midplane. The dashed 
vertical line marks the time at which a large fishbone occurred, with a corresponding 
reduction seen in each of the FI diagnostic signals. 


FIDAsim code calculates weight functions based on the contribution of Monte Carlo 
marker particles, representing reneutralised FI, to the FIDA spectra m- 

With the velocity-space sensitivity of each of the FI diagnostics being determined, 
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Figure 3: Cross-section of the two branches of the D-D fusion reaction as a function 
of the centre-of-mass energy of the reactants. The empirical formula of Bosch & Hale 
(1992) is used to obtain each of the curves [15] . 


diagnostic signals may be modelled based on synthetic FI distributions. These ‘synthetic 
diagnostics’ rely on measured profiles of plasma density and temperature, rotation 
velocity and impurity density. Such modelling forms the basis of section [6} Before 
presenting the experimental results however, we turn our attention to the classification 
of energetic-particle-driven MHD instabilities typically observed in MAST discharges. 

3. MHD activity observed in MAST 

Figure [5] shows a spectrogram of the perturbations to the poloidal magnetic field in 
a typical beam-heated MAST discharge. This is derived from a magnetic pick-up coil 
located on the outboard midplane of the vessel, close to the interior wall. Neutral beam 
injection at E & = 71keV started at 0.05 s in this shot, injecting 2.0 MW of NBI power. 
A second beam at = 61keV was added at 0.18 s, providing an additional 1.5 MW of 
NBI power. 

The period shown in figure [5] exhibits closely spaced bursts of activity close to the 
core TAE frequency for approximately the first 20 ms. The frequency of this activity 
gradually decreases over the following 20 ms, and the bursts start to become more widely 
spaced in time. The low-frequency component of the well-separated bursts between 
0.14 s and 0.18 s appears to span most of the frequency range from the core rotation 
frequency / rot up to the TAE frequency /tae; closer inspection of the signal reveals 
these to be chirping modes , which sweep down rapidly in frequency over approximately 
lms. A clear transition occurs shortly before 0.19 s, when the duration of each burst 
lengthens substantially to approximately 3 ms and the amplitude of the bursts increases 
significantly. These large-amplitude bursts exhibit clear separation into a low-frequency 
fundamental component (n = 1), chirping down from approximately 20 kHz above the 
core rotation frequency to match / rot at the end of the burst, and a second harmonic 
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Figure 4: (a) High-energy part of the FI distribution from shot #29976 at t = 0.2 s, 
R = 1.00 m on the outboard midplane, modelled by Transp. (b) FIDA velocity- 
space sensitivity function (‘weight function’) at A = 661.3 nm for a toroidally-viewing 
chord intersecting the neutral beam at R = 1.00 m. The spectral radiance is obtained 
by multiplying the FI distribution (panel (a)) by the weight function (panel (b)) and 
integrating over velocity space, (c) and (d) are similar to (a) and (b), but for a chord 
intersecting the neutral beam at R = 1.25 m. Note the different ranges of pitch in (c) 
and (d) than in (a) and (b). 


(n = 2) at twice the fundamental frequency in the lab frame. Higher harmonics are also 
observed, albeit at much lower amplitudes. Finally, the bursts undergo a transition to 
a continuous mode which tracks the core rotation frequency from 0.27 s until the end of 
the discharge. 

The evolution of the bursts of MHD activity is seen clearly in a time trace of the 
Mirnov pick-up coil signal. Figure [6] shows the raw data from the coil used to generate 
the spectrogram shown above. There is no clearly defined criterion to establish on 
the basis of the observed behaviour whether a given burst of MHD activity should be 
classed as a TAE, which is a normal mode of the background plasma, or an energetic- 
particle mode (EPM), the existence of which is dependent on the presence of a FI 
population. The classification indicated in figure [6] is based mainly on the frequency of 
the mode. As seen in figure [5j the first set of weakly chirping MHD bursts is clustered 
close to the core TAE frequency, so these modes are identified as TAE. The modes 
in the blue region in figure [6] exhibit much stronger frequency chirping and depart 
significantly from the TAE frequency, suggesting a strong energetic-particle drive which 
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Figure 5: The temporal evolution of the spectrum of poloidal magnetic field 
perturbations on the outboard midplane of the MAST vessel during shot #29207. Black 
and red points show the on-axis toroidal rotation frequency and central frequency of 
the TAE gap respectively. The parameters used to derive these values were obtained 
from charge-exchange measurements of the rotation velocity, Thomson scattering 
measurements of the electron density, and a reconstruction of the magnetic equilibrium 
with the Efit equilibrium solver, /tae plotted here is equal to the sum of the calculated 
core TAE frequency and the toroidal rotation frequency, accounting for the Doppler shift 
of the frequency observed in the lab frame. 


allows them to chirp down into the Alfven continuum panaiza. Features of these 
modes in the Mirnov coil trace include large ‘spikes’ in amplitude lasting for just a few 
wave cycles; asymmetric evolution in time; highly variable spacing between bursts; and 
variable amplitude from one burst to the next. These are tentatively identified as TAE 
with nonlinear coupling to EPM, although they will be referred to simply as ‘chirping 
TAE’ hereafter. In the red region, although some of the features such as the temporal 
asymmetry of the TAE-EPM are still apparent in the first few bursts, the evolution 
of each burst is slower and more regular. The amplitude of successive bursts seems to 
follow a slowly evolving envelope, reaching a maximum at approximately 0.235 s, and 
there is much less inter-event variability. In accordance with the nomenclature used 
for such MHD activity on other devices mmmm, these are called fishbones. A 
characteristic of the fishbones is that the frequency at the end of the burst coincides 
almost exactly with the core toroidal rotation frequency. Finally, the steady-state mode, 
which in MAST tends to evolve directly from the last of a set of fishbones, tracks the 
core rotation frequency and exhibits no frequency chirping and only weak amplitude 
modulation. This behaviour suggests a saturated ideal-MHD mode. Analysis with soft 
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Figure 6: Time trace of the data acquired during shot #29207 by an outboard midplane 
magnetic pick-up coil. The coil is the same as that used to generate the spectrogram 
in figure [5} Overlaid coloured regions show the qualitative classification of the different 
types of MHD activity as it evolves throughout the shot. The classification is based on 
the temporal evolution of the frequency and amplitude of the signal. TAE=toroidal 
Alfven eigenmode; EPM=energetic particle mode; FB=fishbone; LLM=long-lived 
mode, which is the name commonly given to the saturated internal kink mode in MAST. 


X-ray detectors and a toroidal array of Mirnov coils shows the dominant structure of 
this mode to be that of a (ra,n) = (1,1) internal kink [25] - This is often called the 
long-lived mode (LLM) in MAST, due to its tendency to persist from the end of the 
fishbone period until the end of the discharge. 

Figure [7] shows the structure of typical n = 1 TAE and internal kink eigenmodes in 
MAST plasmas, calculated with the linear ideal-MHD code MISHKA-1 [25] - The broad 
structure of the global TAE, with multiple poloidal harmonics peaking at successively 
larger radii, contrasts with the predominantly m = 1 kink mode. The kink mode 
exhibits a ‘double layer’ due to the non-monotonic g-profile, with the internal region 
of negative displacement having q > 1 and the outer region of positive displacement 
having q < 1. Two q = 1 rational surfaces are present, at s ~ 0.2 and s ~ 0.6, which 
form the boundaries of the regions of opposite displacement. In the case of the TAE, 
by contrast, no q = 1 surface exists within the plasma. 

For the purposes of this study, two sets of discharges were selected for the analysis 
of the effects of fishbones and TAE on the FI population. Both sets consist of discharges 
with high NBI power and high electron density, the combination of which provides a 
reasonable signal-to-noise ratio (SNR) in the FI diagnostic signals along with large- 
amplitude bursts of TAE and fishbones separated in time by approximately 5 — 10ms. 
These points are important in allowing the effects of individual bursts of MHD activity 
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Figure 7: a) Structure of the n = 1 TAE calculated with the linear ideal-MHD code 
MISHKA-1 for MAST shot #29210 at t = 0.165 s. The abscissa s represents the square 
root of normalised poloidal flux, and the ordinate represents sv r , where v r is the radial 
velocity of the perturbation in arbitrary units. Integer labels indicate the poloidal mode 
number of the corresponding coloured line, b) As for panel (a), but for the n = 1 
internal kink mode in shot #29976 at t = 0.210 s. 


on the FI population to be resolved by each diagnostic. A further reason for choosing 
these discharges is the fact that the NC was scanned in tangency radius between repeated 
shots, allowing composite NC count-rate profiles to be derived. Global parameters of 
representative discharges from each set are presented in the following section. 

4. Plasma scenarios 

Figure [8] shows time traces of key parameters related to the FI population and MHD 
activity in a representative discharge from each of the two sets analysed in the course 
of this investigation. 

Each set of discharges featured periods with well-separated bursts of MHD activity, 
with the typical evolution of the magnetic spectrogram exhibited by beam-heated 
MAST plasmas in which chirping energetic-particle modes grow in amplitude and 
decrease in frequency as the g-profile evolves; this behaviour was illustrated in section 
|3j Spectrograms from the periods of interest in each of these scenarios are shown in 
figure |9} Shot #29207 et seq. included a period with a single neutral beam injecting 
2.5 MW of power with a primary energy of 71keV, throughout which chirping TAE 
persisted, followed by a period with both neutral beams injecting a total of 3.5 MW of 
power with the primary energy of the second beam set at 61 keV, in which large quasi- 
periodic fishbones were observed to merge into the LLM. The second set of discharges, 
commencing with #29975, exhibited very strong fishbones occurring with semi-regular 
spacing between bursts; the TAE are less distinct, and the early stages of the three 
discharges less similar to each other, than those in the first set, so the analysis of this 
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Figure 8: Global parameters of representative discharges from the two sets selected for 
analysis of the effects of TAE and fishbones. The first set consists of shots #29207 — 
#29210, and the second set includes shots #29975, #29976 and #29980. The panels 
show, from top to bottom: line-integrated electron density measured with a CO 2 laser 
interferometer; total injected NBI power; global neutron rate measured with the fission 
chamber; and RMS amplitude of an outboard-midplane Mirnov coil signal sensitive to 
changes in the poloidal field. 


set focuses on the large fishbones occurring late in the discharge. The total NBI power 
in the second set was 3.0 MW, and the primary beam energies were 70keV and 44keV. 

We now examine the data acquired with the fission chamber, neutron camera 
and FID A diagnostic during each set of discharges. A systematic analysis allows firm 
conclusions to be drawn regarding the effects of chirping modes on the FI population. 

5. Effects of chirping modes on confined fast ions 

5.1. Chirping TAE 

As shown in figure [10[ a significant reduction in the NC count rate and FID A radiance 
is observed to coincide with many TAE bursts in the first set of discharges. The timing 
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Figure 9: Spectrograms from an outboard-midplane Mirnov coil signal measuring the 
poloidal component of the perturbed magnetic field, (a) shows a period during shot 
#29975 in which chirping TAE evolve into large-amplitude fishbones, (b) and (c) show 
two separate periods during shot #29207 in which chirping TAE and large-amplitude 
fishbones respectively dominate the low-frequency magnetic activity. Note the extended 
vertical axis in (b) compared to (a) and (c). In neither of these shots are there any signs 
of high-frequency (compressional) Alfven eigenmodes, which would appear at around 
1 MHz; the absence of high-frequency MHD activity in the presence of low-frequency 
EPM is a characteristic feature of MAST plasmas. 


of these drops in signal relative to the amplitude of the envelope of the Mirnov coil 
signal varies substantially between events. Not all drops in the diagnostic signals are 
associated with prominent chirping modes, neither are all bursts of magnetic activity 
correlated with significant changes in the diagnostic signals. In light of this variability, a 
quantitative approach must be adopted to determine the degree to which the reduction 
in FI confinement indicated by the NC and FIDA signals is correlated with chirping 
TAE. No CFPD data are available for this set of discharges, since the diagnostic had 
not been installed at the time. 

The first stage in identifying the correlation between TAE and FI diagnostics signals 
is to establish criteria for the identification of an ‘event’ in each case. Inspection of time 
traces such as those shown in figure [TO] suggests that strong effects on the diagnostic 
signals are associated with TAE bursts with a large amplitude in the RMS Mirnov coil 
trace. Any increase in magnetic activity which exceeds a threshold RMS amplitude 
is therefore chosen to constitute an ‘event’ in the magnetics signal; this threshold is 
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Time (s) 

Figure 10: a) Traces of NC count rates from the two midplane chords during a period 
in shot #29207 in which chirping TAE were active. The tangency radii of the two 
chords are labelled, and the magenta trace shows the RMS amplitude of an outboard- 
midplane Mirnov coil signal in arbitrary units, b) Traces from two channels of the 
toroidally-viewing FIDA system during the same discharge, with the beam intersection 
radii labelled. Data are averaged over wavelengths 660.7 nm to 661.3 nm, and therefore 
include FIDA emission from passing FI with energy E > 46keV. Error bars on each 
trace represent random uncertainty. 


indicated by the dot-dashed horizontal lines in figure 10 In order to exclude drops in 
the NC and FIDA signals which occur due to statistical noise, or spurious drops in the 
nett FIDA signal caused by rapid, transient spikes in background emission correlated 
with TAE bursts, the condition imposed on these signals is that an ‘event’ constitutes 
a drop in signal which does not recover to the pre-event level within four time steps. 
With the integration time of the FIDA data set at 0.9 ms and that of the NC data set 
at 1.0 ms, this provides sufficient temporal resolution to allow each of the TAE bursts 
to be separated while ensuring that most of the ‘events’ represent true reductions in the 
part of the FI population to which the diagnostics are sensitive. Examples of qualifying 
events in the NC and FIDA signals are shown in figure [TT} Note that a quadratic trend is 
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Figure 11: Drops in NC and FIDA signals which qualify as ‘events’ for the purposes 
of determining the correlation with chirping TAE. The dashed vertical black and blue 
lines mark the onset times of events in the FIDA and NC signals respectively. Dotted 
vertical magenta lines mark the onset times of events in the RMS Mirnov coil signal, 
in which the magenta trace exceeds the threshold shown by the dashed horizontal line. 
The shaded regions which extend from one time step before to two time steps after the 
events in FI diagnostic signals indicate the period within which the magnetics event 
must occur in order to be counted towards successful identification; clearly both of the 
events highlighted here satisfy this criterion. 


removed from each of the time traces prior to analysis; this ensures that significant drops 
in diagnostic signals are not masked by the trend of the signal throughout the analysis 
period, which extends from 0.13s to 0.18s in each of the four discharges analysed in 
this work. 

Proceeding from the strict definition of an ‘event’ in the magnetics and FI diagnostic 
time traces, the criterion chosen to determine whether events coincide is that the onset 
time of the magnetics event must he within a period extending from one time step 
before to two time steps after the start of the drop in FIDA or NC signal. This period 
is indicated by the shaded regions in figure 11 The natural variation between events, 
arising due to the variable levels of drive and damping of the TAE, is thereby taken 
into account. For each FI diagnostic time trace, and for the analysis period during each 
discharge, an ‘identification factor’ may be defined to quantify the correlation between 
events in the diagnostic signals and the occurrence of TAE bursts. This factor is defined 
as 

identified events 


Q ID — 


total TAE events + false positives 


(i) 


The top panel of figure 12 shows the results of the TAE identification analysis for all 
four discharges in this set. In each case, the FI diagnostics are best able to identify the 
effects of TAE bursts in the radial range where the signal is strongest; in the case of 
the NC this corresponds to the region close to the magnetic axis, and in the case of the 
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FIDA diagnostic this corresponds to mid-radius on the outboard side. This observation 
suggests that the limiting factor in the ability of the diagnostics to resolve the effects of 
chirping TAE on FI, at least in this set of discharges, is the SNR. 



0.6 


0.8 


1 


1.2 


Radius (m) 


Figure 12: Top panel: TAE identification factors, defined according to equation [lj for the 
NC and FIDA diagnostics during the period 0.13 s to 0.18 s in shots #29207 — #29210. 
Vertical dashed and dotted lines mark the approximate positions of the magnetic axis 
and q m [ n surfaces respectively. Dashed horizontal lines of various colours mark the 90% 
confidence level for event identification in the corresponding FI diagnostic signal; data 
points have less than a 10% chance of lying above these lines assuming random time 
series. Bottom panel: as above but for fishbones occurring in the period 0.208 s to 0.258 s 
during the same set of discharges. The two innermost FIDA channels are excluded from 
this panel due to the contamination of the spectra by SW beam emission. 


A meaningful assessment of the significance of these identification factors requires 
confidence intervals to be determined. The process adopted here was to randomly 
permute the data points in each time series, and to establish an identification factor Qj D 
based on the new, randomised time series. Repeating this process many times allows 
confidence levels to be quantified as percentiles of all QJ D . If the Qm of the original time 
series lies above, say, the n% confidence level, then the probability of the observed value 
of Qid arising purely by chance is less than (100 —n)%. Dashed, horizontal lines in figure 


level. The advantage of this procedure is that it is insensitive to the distribution of 
data points in the time series as long as the properties of this distribution are constant 
throughout the analysis period. To ensure that this was the case, the time series (with 


11 mark the 90 th percentile of 10,000 values of Q? D , which is to say the 90% confidence 
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least-squares fitted quadratic trends removed) were subjected to the Dickey-Fuller test 
for the presence of unit root m- In all cases, the null hypothesis (that a unit root is 
present in the time series) was rejected at the 1% level, allowing the time series to be 
treated as stationary. It can therefore be said with confidence that chirping TAE cause 
a significant reduction of the confined FI density in these plasmas. 


5.2. Fishbones 


Each of the discharges in the set analysed above also exhibited well-spaced, quasi- 
periodic fishbones as seen in panel (c) of figure [9j The method described in the previous 
subsection was applied to the NC and FID A signals to establish the correlation between 
these signals and the fishbones. In this analysis, the threshold applied to the RMS 
Mirnov coil signal was increased to account for the larger amplitude of the fishbone 
bursts compared to the TAE bursts. Dashed, horizontal lines in each panel of figure [13] 
indicate the chosen threshold. The width of the time window which was searched for a 
burst in Mirnov coil signal corresponding to each drop in NC or FID A signal was also 
widened compared to that used in the TAE analysis. This widening accounts for the 
fact that the fishbones evolve more slowly than the TAE and therefore take longer to 
have a significant effect on the FI population. The window was extended by half a time 
step in each direction, thus covering the period from 1.5 time steps before to 2.5 time 
steps after the commencement of the drop in FI diagnostic signal. Significant drops in 
NC and FID A signals are clearly seen to accompany the fishbone bursts in figure [13] 
The results of the analysis described above are shown in the bottom panel of figure [I2j 
As in the case of the TAE, both diagnostics see significant drops in signal correlated 
with the fishbones. In many of the channels the confidence in the correlation exceeds 
90%. Note that the two innermost channels of the FID A system cannot be included in 
this analysis since the spectra are strongly contaminated with SW beam emission. 

In previous studies of MAST FIDA data [28], it was found that the relative 
change in core FIDA signal due to the fishbones correlated strongly with the maximum 
RMS amplitude of the Mirnov coil signal and only weakly with the amplitude of the 
perturbation. This analysis was repeated for the data from the recent set of discharges 


(#29207 et seq.). Figure 14 shows the scaling of the drops in FIDA signal with the 
amplitude of the fishbones and with the RMS amplitude of the perturbation in the 


Mirnov coil signal. The wavelength used in figure 14 corresponds to a minimum FI 
energy of 46 keV. The data are averaged over three channels with beam intersection radii 
1.03 m, 1.07 m and 1.14 m. A different symbol is used for each shot, but the regression 
analysis takes into account the whole data set. As in the earlier study [28], the drops in 
FIDA signal are correlated much more strongly with the RMS Mirnov coil signal than 
with the fishbone amplitude. A similar linear scaling was found between the drops in 
global neutron rate and the peak Mirnov coil signal in a previous study of MAST data 
m- This suggests that the strength of the interaction between the fishbones and the 
FI is sensitive to the frequency of the perturbation as well as to its amplitude, since 
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a) 



Time (s) 

Figure 13: a) Traces of NC count rates from the two midplane chords during a period 
in shot #29207 in which fishbones were active. The tangency radii of the two chords 
are labelled, and the magenta trace shows the RMS amplitude of an outboard-midplane 
Mirnov coil signal in arbitrary units, b) Traces from two channels of the toroidally- 
viewing FID A system during the same discharge, with the beam intersection radii 
labelled. Data are averaged over wavelengths 660.7 nm to 661.3 nm, and therefore 
include FIDA emission from passing FI with energy E > 46keV. Error bars on each 
trace represent random uncertainty. 


Bq ^ luBq. Here, Bq is the time derivative of the poloidal magnetic field, to which the 
Mirnov coil signal is proportional assuming a fixed mode location, and uj and Bq are 
the frequency and amplitude of the mode. Such a relationship was previously found 
between drops in the global neutron rate and (F?6>)rms during TAE bursts in DIII-D 
m , while no linear dependence on (Bq) rms was seen in the case of fishbone bursts. 
A linear correlation was later found, albeit in a smaller data set, between drops in the 
neutron rate and the amplitude Bq of chirping modes in the fishbone-EPM frequency 
range m- In NSTX meanwhile, a very weak linear scaling was identified between FI 
losses and mode amplitude Bq/B in the case of both EPM and TAE, with TAE having 
a stronger effect for a given mode amplitude m ■ Fishbones were not included in that 
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Fishbone amplitude (io 5 T) Peak RMS Mirnov coil signal (V) 


Figure 14: Relative change in FID A signal plotted against (a) the amplitude of the 
fishbone perturbation and (b) the maximum RMS amplitude of the Mirnov coil signal 
for a set of 23 fishbones in shots #29207 — #29210. Each data point corresponds to 
a single fishbone, and different symbols are used for each discharge. The relative drop 
in signal is defined as (Si — Sf)/Si where Si and Sf are the signal before and after the 
fishbone respectively, The correlation coefficient r, reduced x 2 and cumulative chi-square 
probability P = Pr(X 2 > x 2 ) are shown in each panel. Data are averaged over R = 
1.03 m to 1.14 m at A = 660.7 nm. 


data set, and the correlation with Be was not examined. For the changes in FID A signal 
associated with the chirping TAE discussed in the previous subsection, no meaningful 
correlation was found with either (Bq)r ms or with Bq. The present result therefore raises 
the question of which quantity is the better proxy for the ability of chirping modes to 
redistribute FI, particularly in the presence of mixed low and high-frequency modes. 

Since the tangency radius of the NC was scanned between discharges in this set, 
the correlation analysis presented in figure [14] unfortunately cannot be extended to this 
diagnostic. Furthermore, the drops in global neutron rate due to these fishbones are 
small and difficult to distinguish from noise or from perturbations caused by changes 
in bulk plasma parameters or non-resonant MHD modes. We therefore move on to 
consideration of the possibility that as well as causing real-space or velocity-space 
transport of resonant FI as suggested by the results in this section, the TAE and 
fishbones may also cause losses of FI from the plasma. 


5.3. Fast-ion losses due to resonant MHD instabilities 


Spikes in the D a radiation from the edge of the plasma and from the divertor are 
observed to coincide with many chirping modes. Figures 15 and 16 show examples of 
this behaviour from shot #29207. Bursts in the Mirnov coil trace coincide with a 
significant reduction in FID A and neutron emission from the core of the plasma, as 
well as with pronounced spikes in the signal from a filtered Dq, monitor viewing the 
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E min =55 keV 
E min =58 keV 


E min =61 keV 
E min =64 keV 
E min =67 keV 




Figure 15: Chirping TAE during shot #29207 are observed to cause large drops in core 
NC and FID A signals which in most cases take several milliseconds to recover to the 
pre-event level (top panel). Coincident with the largest of these bursts, spikes in D a 
radiation from the edge of the plasma are observed (middle panel). A large drop in 
FIDA signal and a spike in edge D Q are also observed at around 0.115 s, despite the fact 
that only a small TAE, and no drop in NC signal, is observed at this time. The nature 
of this event is unknown, but the observations suggest a spatially-localised loss of FI. 
The bottom-right panel shows the signal from a background-viewing FIDA channel at 
wavelengths corresponding to minimum detectable FI energies E m [ n close to the beam 
injection energy, which was 71 keV. This signal comprises bremsstrahlung and passive 
FIDA emission, as discussed in the text. The midplane intersection radius of the LOS 
from which these data are obtained is R = 1.33 m. 


edge of the plasma close to the midplane. These spikes in D a light were hypothesised 
to be the result of CX reactions between FI ejected from the plasma by the modes, 
and neutrals close to the plasma boundary. The open divertor in MAST meant that 
deuterium recycled from the plasma at the divertor plates was able to circulate freely 
in the vacuum vessel, forming a dense ‘blanket’ of neutrals at the edge of the plasma. 
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Figure 16: Fishbones during shot #29207 are observed to cause large drops in core NC 
and FID A signals which in most cases take several milliseconds to recover to the pre¬ 
event level (top panel). Coincident with these bursts, spikes in D a radiation from the 
edge of the plasma are observed (bottom panel). The D a spikes in this case exhibit more 
structure, last longer, but have a smaller amplitude than in the case of the TAE. As in 
the previous figure, the bottom right panel shows background emission (bremsstrahlung 
and passive FIDA) at wavelengths corresponding to minimum detectable FI energies 
Emm close to the beam injection energy. 


To test this hypothesis regarding the origin of the D a spikes, the passive FIDA 
emission provides useful information. Passive FIDA is the D a light emitted by 
reneutralised FI which have undergone CX with edge neutrals rather than with beam 
or halo neutrals. By inspecting the time trace of passive FIDA during a time window in 
which edge D a spikes are observed, the mechanism giving rise to these spikes is firmly 
established to be related to the FI. The bottom panels of figures [15] and [16] show time 
traces of passive FIDA from a reference channel (R = 1.33 m) during the period in 
which large-amplitude TAE or fishbones and D a spikes are observed. The spikes in 
passive FIDA emission are clearly seen to be suppressed close to the beam injection 
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energy, which is 71 keV during this shot. This observation demonstrates conclusively 
that the D a spikes are at least partially caused by reneutralised FI rather than by 
thermal neutrals or perturbations to bremsstrahlung radiation. This implies that both 
TAE and fishbones caused losses of FI from these plasmas. 

Analysis and interpretation of the FID A and NC observations reveals that both 
TAE and fishbones cause FI to be transported from the core of the plasma and even 
lost through the plasma boundary. The next step in this analysis is to attempt to model 
the effects of the modes on FI using the ad hoc models available within Transp. 

6. TRANSP modelling of fast-ion redistribution 

6.1. Anomalous transport models 

Two mechanisms are available within Nubeam to allow the effects of TAE and fishbones 
on the FI to be modelled. The redistribution may be modelled as a combination 
of diffusive and convective transport, with user-specified diffusion and advection 
coefficients specified as a function of radial position, time and energy. Without first- 
principles modelling of the underlying physics to guide the choices of these coefficients 
however, only the simplest models are practically applicable. In the modelling described 
in this section, anomalous diffusion was varied as a function of time but was taken to be 
radially uniform throughout the plasma. Convective transport was not applied, since 
diffusive transport was found to reproduce the observations reasonably well. The effect 
of diffusion is to cause radial transport of FI at a rate proportional to the radial gradient 
of the FI density. The resulting flux is 

Tfi = — AmV r n F i, 

where D an is the spatially and temporally-varying anomalous FI diffusivity. The diffusive 
model was applied successfully to model the global neutron rate and NC count rate 
profiles in MAST discharges with on- and off-axis NBI [29] . 

The second FI transport mechanism available within Nubeam is an ad hoc fishbone 
loss model. This model takes as input parameters: the time of the first fishbone and 
the time at which the model is switched off; the duration of the fishbone burst and the 
characteristic loss time of the FI while the fishbone is active; the time between fishbone 
bursts; and the minimum and maximum values of energy and pitch of the affected FI. 
This model was successfully used to reproduce global and local neutron emission from 
MAST plasmas as described by Klimek et al. [Mj. In the present work however, FIDA 
observations are used as an additional constraint on the parameters of the model. 

In the following subsections the results of Transp modelling using anomalous 
diffusion and the fishbone model are presented. 

6.2. Modelling with anomalous diffusion 

The series of discharges #29207 et seq ., with chirping TAE and quasi-periodic fishbones, 
provided an opportunity to test the diffusive Nubeam modelling in the presence of 
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resonant MHD activity. Time slices averaged over 3 ms, from 0.162 s to 0.165 s during 
the chirping TAE phase and from 0.237 s to 0.240 s during the fishbone phase, were 
used for a comparison between Nubeam modelling and FI diagnostic signals under the 
assumption of spatially-uniform anomalous diffusion. The results of this comparison are 
shown in figure [17] 

Clearly the introduction of anomalous diffusion is sufficient in this case to model 
the effects of both chirping TAE and fishbones on the FI distribution, at least to the 
extent that the distribution may be constrained by the available diagnostics. It is worth 
noting that at the two times shown in figure [17] the measured profiles are averaged over 
periods during which a significant MHD perturbation occurred; profiles during selected 
MHD-quiescent windows shortly before a TAE or fishbone burst exhibit increased signal 
levels, although the signal still does not reach the level modelled in the absence of 
anomalous FI diffusion. The coarse assumption of spatially-uniform diffusivity, held 
constant over a long period, cannot capture the dynamics of individual bursts of MHD 
and the interaction of these bursts with the FI. 


6.3. Inclusion of the fishbone model 


Figure [18] shows that spatially-uniform anomalous diffusion models well the early part of 
shot #29975 et seq. during which chirping TAE were active. The forward-modelled NC 
and FIDA profiles show good agreement with the measured profiles at t = 0.160 s, when 
a diffusivity of D an = 3.0 m 2 s -1 allows the Transp modelling to match the measured 
global neutron rate to within 10%. 

Later in each of the discharges from this set, a series of large fishbones caused 
transient FI redistribution which the simple diffusive transport model is unable to 
capture. This redistribution resulted in drops of ~ 10% in the global neutron rate, 


as seen in figure 19 As mentioned earlier, Nubeam contains a fishbone model which 
allows FI to be expelled from the plasma within a specified part of velocity space, with 
a chosen periodicity and characteristic loss time. This model was previously used to 
model changes in global and local neutron emission in another series of MAST discharges 
[55] . but FIDA data were unavailable for those shots. In the series of discharges studied 
here, the fishbone model was initially used with settings similar to those applied in 
the previous study, cutting out all FI with E > 50keV and with —0.5 < p < 0.5. It 
was found however that despite the relatively successful match to the global neutron 
rate and the changes in NC count rate profiles achieved using this model, the changes in 
FIDA signal were not reproduced. This may be understood by referring to the geometric 
weight functions of the FIDA system discussed in section [2j The toroidal FIDA views 
are insensitive to ions with such a small ratio of parallel to perpendicular velocity, since 
the LOS are almost tangential to the magnetic field at the point of beam intersection in 
the core of the plasma. By contrast, the large fishbones in these shots are observed to 
cause significant drops in toroidal FIDA signal at all radii from the magnetic axis to the 
plasma boundary. This observation informed the choice of energy and pitch boundaries 
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Figure 17: Comparisons between measured and Nube AM-modelled FI diagnostic time 
traces and radial profiles in the presence of chirping TAE and fishbones. Panel (a) shows 
the global neutron rate measured in shot #29210 (black) compared to that modelled 
with (red) and without (blue) spatially-uniform anomalous FI diffusion in Transp. The 
values of anomalous diffusivity used in these simulations are given in the labels. The 
forward-modelled NC and FID A profiles are compared to the measured profiles in panels 
(b) and (c) for t = 0.165 s and in panels (d) and (e) for t = 0.240 s. Contamination of the 
core FIDA signal with SW beam emission at the later time means that the two channels 
between 0.9 m and 1.0 m cannot be included, and the signal in the two channels between 
1.0 m and 1.1m is artificially elevated. Vertical dashed and dotted lines in panels (b) 
to (e) indicate the approximate radial positions of the magnetic axis and minimum-g 
surfaces. A = [660.9 — 661.9] nm for FIDA data. 


in the fishbone model within which the FI were removed. An iterative process was 
adopted by which the simulation parameters were at first coarsely modified to converge 
on the global neutron rate, then modified more precisely to converge on the relative 
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Figure 18: Profiles of measured and modelled NC and FID A signals averaged over the 
period 0.157 s to 0.160 s during shot #29976 (with the measured NC profile including 
data from similar discharges #29975 and #29980). The modelled profile is based on a 
simulation with an anomalous diffusivity of 3.0 m 2 s -1 at this time during the discharge. 
Vertical lines mark the approximate radial positions of the magnetic axis and minimum-g 
surface. 


changes in NC and FIDA radial profiles. The neutron rate comparisons from the two 
most successful simulation runs are shown in figure [19] 

Figure 20 shows that Transp run 070, in which high-energy FI with 0 < p < 0.7 
are removed by the fishbones, produces a good match to the relative size of the observed 
changes in NC signal across the plasma radius as well as to the shape of the NC profile. 
Run 069 produces a similar match to the profile shape, but significantly underestimates 
the size of the drops in signal due to the fishbones. 

The same pair of Transp runs was used to generate synthetic FIDA profiles. 
Although the absolute match between the modelled and measured signal levels was 
slightly better for run 070, the simulated fishbones with these settings failed to cause 
the observed reduction in core FIDA signal. Figure |2T] shows that run 069, by contrast, 
matches the relative change in signal fairly well and also matches the profile shape well, 
at least at mid-high energies (top panels), even though it underestimates the magnitude 
of the signal. The comparison between the two runs in terms of their match to the 
relative change in FIDA signal is quantified in figure 22 Clearly the removal of high- 
energy passing ions reproduces the FIDA observations. The resonant interaction of 
passing FI with fishbones finds a precedent in results from the PBX tokamak [33], and 
was explained theoretically by Betti and Freidberg m- 

The fact that run 069 systematically slightly overestimates the magnitude of the 
changes in FIDA signal across the plasma radius, while run 070 underestimates these 
changes in the core but matches them well at the edge, suggests that the real effect of 
the fishbones is best represented by a combination of these models. It is reasonable to 
suppose that the properties which allow FI to resonate with the mode, and hence to be 
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Figure 19: A comparison between the measured (black trace, with yellow shaded 
region indicating statistical uncertainty) and modelled (blue and magenta traces) global 
neutron rates for shot #29976 with the fishbone model applied within Nubeam. The 
fishbones were specified to occur from 0.201s to 0.230 s with a period of 12.8 ms and a 
duration of 4 ms. The characteristic loss time of the affected FI was 3 ms in run 069 
and 1ms in run 070. In run 069, FI with 50keV < E < 75keV and 0.69 < p < 0.93 
(evaluated at either of the two previous midplane crossings) were removed by the 
fishbones; these values of pitch correspond to co-passing FI. In run 070, these ranges 
were 60keV < E < 70keV and 0 < p < 0.7, corresponding to trapped and co-passing 
FI. In addition to the fishbone model, a small anomalous diffusivity of 0.15m 2 s _1 was 
applied to FI with E < 45keV during the fishbone phase in both simulations; D an 
decreased linearly with energy from 0.150 m 2 s _1 to 0.105 m 2 s _1 between 45keV and 
75 keV. Grey and red shaded rectangles indicate the periods, before and after the first 
two fishbones, over which synthetic FI distributions and measured NC and FID A signals 
were averaged to produce the subsequent figures in this section. 


strongly redistributed, are a function of radial position in the plasma. This dependence 
cannot be captured with the simple ad hoc model available in Nubeam. More realistic 
simulations must await the inclusion of a first-principles model of resonant transport in 
the global transport simulations. 

In the set of shots commencing with #29975, the radial position of the CFPD was 
scanned between discharges. This allowed the effects of the large fishbones on the CFPD 
count rate to be determined as a function of midplane intersection radius; the data from 
each channel are shown in figure [23} It is apparent that the strongest effect of the 


fishbones on the count rate is observed on trajectories intersecting the midplane close to 
the magnetic axis. Channels with midplane intersection radii from 0.86 m to 0.96 m see a 
strong depletion in the count rate coinciding with each of the large fishbones, while those 
at 0.79 m and 0.81m see only a weak effect. Without a reliable synthetic diagnostic, 
the fishbone modelling carried out with Transp/Nubeam cannot be validated against 
these measurements. Nonetheless, the observation of fishbone-induced drops in signal in 
all of the traces in figure [23] is broadly consistent with the NC and FIDA measurements, 
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Figure 20: Profiles of NC count rate before (black points) and after (red points) the 
first two fishbones in shots #29975, #29976 and #29980. Panels (a) and (b) compare 
the measured profiles to those forward-modelled from the neutron emission profiles 
generated in Transp run 069 for shot #29976 (dashed lines), while panels (c) and 
(d) show a comparison with Transp run 070 for the same shot. Vertical dashed and 
dotted lines mark the radial positions of the magnetic axis and minimum-# surface. 


which show a significant reduction in signal across the core of the plasma. 

7. Summary and conclusions 

The results presented in this paper demonstrate the advantages of employing several 
complementary diagnostics to study the behaviour of the FI population. Extending 
previous studies on MAST [HI [23 123 ES], a systematic analysis of NC and FID A data 
has allowed the deleterious effects of chirping TAE and fishbones on the confined FI 
population to be identified with confidence. The velocity-space sensitivity of the FIDA 
measurements moreover permits losses of FI to be identified unambiguously, and these 
losses are seen to be strongly correlated with the MHD activity. 

Radial profiles of neutron and FIDA emission have been used to constrain transport 
modelling and to establish the effects of fishbones on the FI distribution. A simultaneous 
comparison of the modelling results to measurements made with the FC, NC and FIDA 
diagnostic has allowed constraints to be placed on the effects of fishbones on the velocity- 
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Figure 21: Profiles of FID A radiance before (black points) and after (red points) the first 
two fishbones in shot #29976. Panels (a) and (b) compare the measured profiles to those 
forward-modelled with FIDAsim from the FI distributions generated in Transp run 
069, in a wavelength range corresponding to a minimum FI energy of 48keV (dashed 
lines). Panels (c) and (d) show a comparison at higher wavelengths, with E m [ n = 57keV. 
Vertical dashed and dotted lines mark the radial positions of the magnetic axis and 
minimum -q surface. 


space distribution of FI; these modes are seen to affect high-energy, passing particles 
in the core of the plasma. Data from the prototype CFPD support the conclusion that 
fishbones have a strong impact on core FI confinement. 

Resonant energetic-particle-driven MHD instabilities observed in MAST, as 
discussed in section [3j include chirping TAE and fishbones. In this work, observations 
of changes in the FI diagnostic signals correlated with each of these instabilities were 
presented and the interpretation of these observations was discussed. The results of 
Transp modelling in which anomalous FI transport or loss processes were invoked to 
try to match the observations were presented and compared with the measurements. 
Four main conclusions are supported by the results presented here: 

• Drops in both FID A and NC signals are correlated with chirping TAE and fishbones 
at a statistically significant level, indicating a reduction of the confined FI density. 

• Chirping TAE and fishbones cause enhanced losses of FI from the plasma. 

• The profiles of NC and FIDA signals averaged over individual bursts of MHD 
activity are well modelled by applying anomalous diffusion to the FI in Nubeam. 
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Figure 22: Profiles of the relative reduction in FIDA radiance, defined as (Si — Sf)/Si , 
where Si and Sf are the pre and post-fishbone signals, caused by two fishbones in 
shot #29976. The measured relative changes are compared to those derived from the 
FIDAsim profiles based on Transp run 069 (blue dashed lines) and 070 (green dashed 
lines). Wavelengths are the same as those in the corresponding panels of figure 
panels (a) and (c) correspond to the first of the modelled fishbones, while (b) and (d) 
correspond to the second fishbone. 
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• Fishbones strongly affect the high-energy, passing FI population in a manner which 
may be reproduced with some success using the fishbone model in Nubeam. 

The additional flexibility offered by the NBI system on the forthcoming MAST- 
Upgrade device, which will have the capability to inject beams simultaneously on and 
off-axis, will allow a detailed investigation into the relationship between beam deposition 
profiles, MHD mode activity and FI transport [55] . 
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Figure 23: Time traces of data from two channels of the CFPD during the large 
fishbones in shot #29975 et seq.. Each column includes data from a single shot, with 
dotted vertical lines marking the time at which the amplitude of each of the three 
largest fishbones reached a maximum. The label in each panel indicates the midplane 
intersection radius of the centre of the bundle of fusion proton trajectories which reaches 
the detector. The magnetic axis is located at approximately R = 0.97 m. 
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